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NONLINEAR STRUCTURE FORMATION, BACKREACTION 
AND WEAK GRAVITATIONAL FIELDS 

Aseem Paranjape^ 

Abstract. There is an ongoing debate in the hterature concerning the 
effects of averaging out inhomogeneities ( "backreaction" ) in cosmology. 
In particular, some simple models of structure formation studied in the 
literature seem to indicate that the backreaction can play a significant 
role at late times, and it has also been suggested that the standard 
perturbed FLRW framework is no longer a good approximation dur- 
ing structure formation, when the density contrast becomes nonlinear. 
In this work we use Zalaletdinov's covariant averaging scheme (macro- 
scopic gravity or MG) to show that as long as the metric of the Uni- 
verse can be described by the perturbed FLRW form, the corrections 
due to averaging remain negligibly small. Further, using a fully rela- 
tivistic and reasonably generic model of pressureless spherical collapse, 
we show that as long as matter velocities remain small (which is true in 
our model), the perturbed FLRW form of the metric can be explicitly 
recovered. Together, these results imply that the backreaction remains 
small even during nonlinear structure formation, and we confirm this 
within the toy model with a numerical calculation. 

1 Introduction 

' The application of general relativity (GR) to the large length scales relevant in 
. cosmology, necessarily requires an averaging operation to be performed on the 
■ Einstein equations (Ellis 1984). The nonlinearity of GR then implies that such an 
averaging will modify the Einstein equations on these large scales. Symbolically, 
this happens since E[{g)] ^ {E[g\), where g is the metric and E the Einstein 
tensor. 

In recent times, attention has been focused on two promising candidates for 
a consistent nonperturbative averaging procedure, namely the spatial averaging of 
scalars (Buchert I2000p and the covariant approach known as macroscopic grav- 
ity (MGHZalaletdinov 119921 11993p . The magnitude of the corrections has been 
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debated in the literature (see e.g. Kolb et al. 120061 Ishibashi & Wald [20061 Rasa- 
nen 120061 Hirata & Seliak [^UD^ . Broadly speaking, due to the structure of the 
Riemann tensor, the modification terms C in any averaging approach, will have a 
symbolic structure given by 

C~(f2)-(f)2, (1.1) 

where T denotes the Christoffel connection, and the tilde represents any processing 
of the Christoffels required by the averaging operation. 

In this talk we will work within the MG framework of Zalalctdinov, since 
this framework allows us to consistently define an averaged metric, which we will 
assume to be of the flat FLRW form 

c^SpLRw = -d'r'^ + a{Tfdx^ . (1.2) 

Specifically, we will work in the spatial averaging limit of MG, which is relevant 
to cosmology (Paranjape & Singh 2007J. We will also assume, consistently with 
observations of the cosmic microwave background (CMB), that the early universe 
was well described by a metric which is a perturbation around the FLRW form, 
given by 

ds^ = -(1 + 2ip)dT^ + a(r)2(l - 2il))dx^ . (1.3) 

The main argument of this talk is that although technically possible, in the real 
world it is extremely unlikely that backreaction significantly influences the average 
cosmological expansion. In section 2 we will detail the calculation of the backre- 
action during the linear regime of perturbation theory, and mention some lessons 
that one can learn from this calculation. In section 3 we will describe a toy model 
of fully relativistic and nonlinear structure formation, and demonstrate that the 
metric for this model can be explicitly brought to the perturbed FLRW form (|1.3p , 
even during the nonlinear phase of the evolution. A calculation of the backreac- 
tion along similar lines as in the linear case then shows that the backreaction in 
this toy model remains small even at late times, thus supporting our argument. 

2 Building the argument 

2. 1 The linear regime 

The expression in (jl.ip can be used to obtain a simple estimate for the backre- 
action, assuming a perturbed FLRW metric of the form (|1.3p with ip — ip in the 
absence of anisotropic stresses. We then have C ^ a~'^{Vip ■ Vi^). Assuming a 
two component flat background consisting of cold dark matter (CDM) and radi- 
ation (known as standard CDM or sCDM), and taking the averaging to be an 
ensemble average over the initial conditions, it is not difficult to show that in the 
matter dominated era, this estimate leads to C ~ [10^^i?Q/a(T)^j3. This indicates 



^ The factor 10 * arises from a product of the normaUsation of the initial power spectrum 
A ~ 10-3, and the factor (keq/Ho)'^ ~ 10= which arises in the transfer function integral, where 
ksq = aeqH{aeq) is the wavenumber corresponding to the radiation-matter equality scale. 
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that at least for epochs around the last scattering epoch, the backreaction due to 
averaging was negligible. 

The real situation is somewhat more complex than this simple calculation 
indicates. On the one hand, the time evolution of a(r) is needed in order to 
solve the equations satisfied by the perturbations, as we effectively did above by 
assuming a form for a(T). On the other hand, the evolution of the perturbations 
is needed to compute the correction terms C. Until these corrections are known, 
the evolution of the scale factor cannot be determined; and until we know this 
evolution, we cannot solve for the perturbations. 

To break this circle, we will adopt an iterative procedure. We first compute 
a "zeroth iteration" estimate for the backreaction, by assuming a fixed standard 
background a^°^ such as sCDM, evolve the perturbations and compute the time 
dependence of the objects C, denoted C'"^. Now, using these known functions of 
time, we form a new estimate for the background a*^^-* using the modified equations, 
and hence calculate the "first iteration" estimate C^'^\ This process can then be 
repeated, and is expected to converge as long as perturbation theory in the metric 
remains a valid approximation. 

With these ideas in mind, we can go ahead and compute the "zeroth iteration" 
estimate for the backreaction in the Zalaletdinov framework. This is given by the 
following equations (Paranjape i2008|l 



1 da 
a dr 
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(2.1) 
(2.2) 

(2.3) 

(2.4) 

(2.5) 

(2.6) 
with 

H = a' /a, $fc is the Fourier space transfer function defined by (p{k, rj) — (p^^^k{il)i 
and Pipi (k) is the initial power spectrum of ip. The results of a numerical calculation 
are shown in figure [H where all functions are normalised by the Hubble parameter 
iJ^(a), and confirm that this zeroth iteration estimate in fact gives a negligible 
contribution. Further, were we to carry out the next iteration, we would essentially 
obtain no difference between the zeroth and first iteration scale factors upto the 
accuracy of the calculation, and hence this iteration has effectively converged at 
the first step itself. 



Here, the prime denotes a derivative wrt. conformal time (' = 9^ = adr 
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Fig. 1. The backreaction for the sCDM model, normalised by -ff^(a). 5*^', p'^' and 5'^' 
are negative definite and their magnitudes have been plotted. The vertical line marks 
the epoch of matter radiation equality a = aeq. 



2.2 Lessons from linear theory 

We see that the dominant contribution to the backreaction at late times, is due 
to a curvature-like term ^ a^^, as expected from our simple estimate above. In 
order to obtain a correction which grows faster than this, we need a nonstandard 
evolution of the metric potential which can only happen if the scale factor 
evolves very differently from the sCDM model, which in turn would require a 
significant contribution from the backreaction. The same circle of dependencies as 
before, now implies that as long as the metric is perturbed FLRW, the backreaction 
appears to be dynamically suppressed. 

Secondly, as figures [5] and [3] show, scales which are approaching nonlinearity, 
do not contribute significantly to the backreaction, which is a consequence of 
the suppression of small scale power by the transfer function. We will return to 
this point when discussing the backreaction during epochs of nonlinear structure 
formation. 



3 The nonlinear regime 

Let us now ask whether one can make meaningful statements concerning the 
backreaction during epochs of nonlinear structure formation, when matter den- 
sity contrasts become very large and perturbation theory in the matter variables 
has broken down. We begin by considering some order of magnitude estimates. 

3.1 Dimensional arguments, and why they fail 

Let us start with the assumption that although the matter perturbations are large, 
one can still expand the metric as a perturbation around FLRW. We are looking for 
either self-consistent solutions using this assumption, or any indication that this 
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Fig. 2. The dimensionless integrand of .S''^', namely the function {k/Ho)'^^k, at three 
sample values of the scale factor. The function dies down rapidly for large k, with the 
value at some k being progressively smaller with increasing scale factor. The declining 
behaviour of the curves for a = Ueq and a = 200aeq extrapolates to large k. 
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Fig. 3. The dimensionless CDM density contrast. Together with figure[3]this shows that 
nonlinear scales do not impact the backreaction integrals significantly. 

assumption is not valid. Given that the metric has the form l|1.3p (and further 
assuming (p = ijj as before), the relevant gravitational equation at late times and 
at length scales small comparable to H~^, is the Poisson equation given by 

^VV = 47rGp(5, (3.1) 

where 6 = {p{t, x)/ p{t) — 1) is the density contrast of CDlvd. As before, we can 
estimate the dominant backreaction component to be C ' Vv'). 



^We are only worried about the small, sub-Hubble scales, since larger scales are well described 
by linear theory where we know the form of the backreaction. 



6 



Title : will be set by the publisher 



Now, for an over/under-density of physical size R, treating a ~ i? ^ and 
Gp on dimensional grounds, we have 

\v\-{HRr\S\ . (3.2) 

For voids, we can set S ~ —1, and then C ^ H^{HR)^ <^ H^, since we have 
assumed HR <C iH . This shows that sub-Hubble underdense voids are expected to 
give a negligible backreaction. For overdense regions we need to be more careful, 
since here S can grow very large. In a typical spherical collapse scenario, the 
following relations hold, 

i? - (1 - cosM)r ; H-^ - (Gp)"^/^ - t - Hf^\u - sinu) , (3.3) 

^ -WR- ^ (1 - C0SU)3 ' ^ ^ ^P'P^ ^ (1 - C0S.)3 ' (3-4) 



which lead to 



(1 — cosu) 



1 



(3.5) 



It would therefore appear that at late enough times, the perturbative expansion 
in the metric breaks down with \ip\ ~ 1, and the backreaction grows large \C\ ~ 1. 
However, the crucial question one needs to answer is the following : Is this situation 
actually realised in the universe, or are we simply taking these models too far? 
We make the claim that perturbation theory in the metric does not break down at 
late times, since observed peculiar velocities remain small. The spherical collapse 
model is not a good approximation when model peculiar velocities in the collapsing 
phase grow large. To support this claim, we will work with an exact toy model of 
spherical collapse. 



3.2 Calculations in an exact model 

The model we consider was presented by Paranjape & Singh (|2008ap . and can be 
summarized as follows. The matter content of the model is spherically symmetric 
pressureless "dust" , and hence the relevant exact solution is the Lemaitre-Tolman- 
Bondi (LTB) metric given by 

ds^ = -dt^ + J^^^^ + R^dn\ (3.6) 
1 — k[r)r'^ 

Here t is the proper time measured by observers with fixed coordinate r, which 
is comoving with the dust. R{t,r) is the physical area radius of the dust shell 
labelled by r, and satisfies the equation R^ — 2GM{r)/R — k{r)r'^ . Here M{r) is 
the mass contained inside each comoving shell, and a dot denotes a derivative with 
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Fig. 4. Density contrast in the LTB toy model, at r = 8.35 Mpc. to = 2/3Ho ~ 9 Gyr. 



respect to the proper time t. The energy density of dust measured by an observer 
comoving with it satisfies the equation p(t,r) — M' {r) / Ai:E? R' , where the prime 
now denotes a derivative with respect to the LTB radius r. 

Initial conditions are set at a scale factor value of — 10^'^, and are chosen 
such that the initial situation describes an FLRW expansion with a perturbative 
central overdensity out to radius r = , surrounded by a perturbative underden- 
sity out to radius r = r„ , with appropriately chosen values for the various param- 
eters in the model (see Table 1 of Paranjape & Singh [^08ap . Figure S] shows the 
evolution of the overdensity contrast in the central region. Clearly, at late times 
the situation is completely nonlinear. Nevertheless, it can be shown that a coordi- 
nate transformation in this model can bring its metric to the form (|1.3p . provided 
one has \av\ <C 1 where v = dt{R/a). Physically v is the "comoving" peculiar 
velocity. The metric potentials (which are actually equal at the leading order, see 
Van Acoleven [2008]) have the expressions Lp = + {l/2){av)'^, ip = £_°H + ^, 
where ^ and are obtained by integrating ^' = (l/2)(fc(r)r^ + {avY){R' / R) and 
^0' = avR' . A numerical calculation shows that av and hence the metric potentials 
do in fact remain small for the entire evolution, for this model. Further, the infall 
peculiar velocity can only become large if the true infall velocity R is large, in 
which case the specific background chosen to define the peculiar velocity, becomes 
irrelevant (since HR <C 1). Hence, the fact that relativistic infall velocities are not 
observed in real clusters etc., leads us to expect very generally that the perturbed 
FLRW form for the metric should in fact be recoverable even at late times. 

Finally, figure [5] shows the dominant contribution to the backreaction in the 
toy model (Paranjape & Singh [2008b[) . There is a significant departure from a 
curvature-like behaviour, due to evolution of the metric potentials. More impor- 
tantly, the maximum value of the backreaction here is 10'^ H^, as opposed to 
^ 10"'* -ff^ as seen in the linear theory. This can be understood by noting that 
the inhomogeneity of our toy model is only on relatively small, nonlinear scales 
(< 20ft,"*Mpc), and the value of the backreaction is therefore consistent with our 
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Fig. 5. The evolution of \S''''^\/6H^ . Also shown is a hypothetical curvature-like correc- 
tion, evolving like ~ a~^. 

earlier observation that nonlinear scales contribute negligibly to the total backre- 
action. To conclude, we have seen that as long as the metric has the perturbed 
FLRW form, the backreaction remains small. Further, there are strong reasons to 
expect that the metric remains a perturbation around FLRW even at late times 
during nonlinear structure formation, a claim that is supported by our toy model 
calculation. It should be possible to test this claim in W-body simulations as well. 
It appears therefore, that backreaction cannot explain the observed acceleration 
of the universe. 

I thank the organisers, especially Thomas Buchcrt, for their kind hospitality 
and for holding a very lively conference in the beautiful city of Lyon. 
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